Calculation of the One- and Two-Loop Lamb Shift for Arbitrary Excited Hydrogenic States 
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General expressions for quantum electrodynamic corrections to the one-loop self-energy [oforderQ(Za)'^] 
and for the two-loop Lamb shift [of order {ZaY'] are derived. The latter includes all diagrams with closed 
fermion loops. The general results are valid for arbitrary excited non-S hydrogenic states and for the normalized 
Lamb shift difference of S states, defined as A„ = n"^ AE(nS) — AE{1S). We present numerical results for 
one-loop and two-loop corrections for excited S, P and D states. In particular, the normalized Lamb shift 
difference of S states is calculated with an uncertainty of order 0.1 kHz. 
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The theory of quantum electrodynamics, when applied to 
the hydrogen atom and combined with accurate measure- 
ments (HQ], leads to the most accurately determined phys- 
ical constants today |3] and to accurate predictions for tran- 
sition frequencies. Of crucial importance are higher-order 
corrections to the bound-state energies, which involve both 
purely relativistic atomic -physics effects and are mixed with 
the quantum electrodynamic (QED) corrections. In general, 
this leads to a double expansion for the energy shifts, both in 
terms of the QED coupling a (the fine- structure constant) and 
the nuclear charge number Z. 

As is well known, the leading one-loop energy shifts (due to 
self-energy and vacuum polarization) in hydrogenlike systems 
are of order a {Z a)* in units of the electron mass. Analytic 
calculations for higher excited states in the order a{Z a)^ are 
extremely demanding. For non-S" states, the a (Z (xf cor- 
rections have been obtained recently |4]. However, excited 
S states are very important for spectroscopy, and the corre- 
sponding gap in our knowledge is filled in the current Letter 
(see Table U. Regarding the two-loop correction, complete 
results for the o? iZ a)^ effect were obtained in 1970 (see 
Ref. 1 5]). Here, we derive general expressions which allow 
the determination of the entire two-loop o? {Z oif' correction, 
for all non-5 hydrogenic states and the normalized difference 
A„ = r? AE{nS) — AE{1S), including the nonlogarithmic 
term. Together with other available analytic |6, 7] and nu- 
merical calculations for the IS state |8], our results allow for 
a much improved understanding of the higher-order two-loop 
corrections for a general excited hydrogenic states, and pave 
the way for an improved determination of fundamental con- 
stants from hydrogen spectroscopy. 

The one-loop bound-state self-energy, for the states under 
investigation here, can be written as 



6(^^E = {A, 



{Zaf [A6iln[(Za)-2]+^g„]}, 



where the indices of the coefficients indicate the power of Z a 
and the power of the logarithm, respectively. We work in 
D = 4—2 e spacetime dimensions, and the dimension of space 
is d = 3 — 2 e. Units are chosen so that ft — c = ~ 1, and 
the electron mass is unity. A nonrelativistic, "Bethe-style" 0] 



calculation of the contribution due to ultrasoft photons, in the 
dipole approximation, leads to a dimensionally regularized 
energy shift E^q, 
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where In fco = 2(Za)i 
is the Bethe logarithm, and 6'^{r) = V'^V/{Ait) is a d- 
dimensional Dirac delta function obtained via the action 
of the Laplacian on the d-dimensional Coulomb potential 
V{r) = -Zar^-'* [r(| - l) 7ri-''/2j^ All matrix ele- 
ments (•) are to be evaluated with regard to the reference 
state, as given by a nonrelativistic (Schrodinger-Pauli) wave 
function, and the summation convention is used throughout 
this Letter 

TABLE I: Values of the nonlogarithmic self-energy correction Aqq 
("relativistic Bethe logarithm") for higher excited S states. 
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AeoinS) 
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A(io{nS) 


1 


-30.924 149 46(1) 


5 


-31.455 393(1) 


2 


-31.840 465 09(1) 


6 


-31.375 130(1) 


3 


-31.702 501(1) 


7 


-31.313224(1) 


4 


-31.561922(1) 


8 


-31.264 257(1) 
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EMM 

, we now consider corrections due to 
the relativistic Hamiltonian, the quadrupole term and the rel- 
ativistic and retardation corrections to the current. The rela- 
tivistic coiTection to the Hamiltonian is 

iJfl = + ^ Z a 5'^{r) + -- a'^ V'V p> . (2) 



Here, cr*-' = [cr', cr^]. The resulting, dimensionally regu- 
larized, correction to the Bethe logarithm is 
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where L{Za) = In [i(Za)-2] , and G = 1/{E ~ H)' is the 
reduced Green function; j3i is a generaHzed Bethe logarithm. 
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We temporarily restore the reference state in the notation of 
the matrix element, and the sums over n and m include both 
the discrete as well as the continuous part of the spectrum. 
The argument of the logarithm in /3i is ln[|iJ — E\/ [Z a)^], 
not ln[2|i/ — E\/{Z a)^] as in In Uq, and this fact is important 
for the precise definition of and of all other generalized 
Bethe logarithms in the following. 

In the dimensional scheme, the quadrupole correction El2, 
which was denoted as F^q in former work , is found 

to be expressible as E]^2 — 1^2 + where 
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and !F2 contains the generalized Bethe logarithm /32, 
a{Zaf[}2 
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Here, n is a three-dimensional unit vector, and we integrate 
over the entire solid angle fi^;. Throughout this Letter, and 



V'' are understood to exclusively act on the quantity immedi- 



ately following the operator, i.e. (V'^V p'^j = ^(V^V^)p^ 
V'^VGH]^^ = ({^^V)GHr^ etc. 

The correction E]^^ to the transition current reads Ej^^^ = 
I^s + ^3, where = a{ZaY P^/t^h^ contains the general 
ized Bethe logarithm (3^, and 
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Here, = p'p^ + i cr'-' V^V, and V = 9/9r* denotes the 
derivative with respect to the ith Cartesian coordinate. The di- 
vergences (in e) in the corrections to the Bethe logarithm are 
compensated by high-energy virtual photons, which in non- 
relativistic QED (NRQED) are given by effective operators. 
From a generalized Dirac equation (see Chap. 7 of Ref. 1 121 ), 
one easily obtains the effective one-loop potential 
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which in leading order gives rise to the correction ^(5'^^)y). 
This correction is a contribution to the middle-energy part 
Em, which originates from from high-energy virtual photons, 
with electron momenta of order Za. The corrections of rela- 
tive order {Z aY to ((J^^-'y) involve relativistic corrections to 
the wave function and to the operators, and a two-Coulomb- 
vertex scattering amplitude. The sum is 

EM = i^^^^^y') + 2(b^^^VGHj^ (6) 

+ ?(l^-4fc)^^^^ + ^^-'>^^^^^^^ 

-3l^({^"'^^^ + 2-^^'^'^^1> 

TT V240 me)^ ' t: V48 3ey \^ ' / 

The complete one-loop result S'^^^'E = Elo + Eli + El2 + 
Ehi + Em reads 
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The matrix elements in this result can be evaluated using stan- 
dard techniques. In terms of the notation of Ref. |4], we have 
L = X]^=i A- O^i" general result Q, evaluated for hydro- 
genic states, reproduces the known logarithmic term Agi, and 
is consistent with all formulas reported for the nonlogarithmic 
term in Eqs. (10) and (12) of Ref. 0]. The evaluation of C, 
is a demanding numerical calculation, and numerical values 
for non-5 states have been presented in Table I of Ref. |4]. 
Taking advantage of the result 1 10] for 15 and the validity of 



Eq. Q for the nS-lS difference, we can now proceed to in- 
dicate results for the nonlogarithmic term Aqq for nS states, 
an evaluation made possible by our generalized NRQED ap- 
proach (see Table|l|i. 

A generaUzation of our NRQED approach leads to the 
following general result for the a"^ (Za)^-term of the com- 
plete two-loop Lamb shift (including all diagrams with closed 
fermion loops). 
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2179 9 5 9 

197 3 1 3 " 
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C(2) ln(2) 
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32 
59 



C(2) 
C(2)- 
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' 64 
17 
32 



C(3) 
C(3) 



87697 9 , 2167, 

+ -C(2) ln(2) - — -C(2) 



345600 10 
16841 1 



207360 



C(2) ln(2) 



9600 
223 
2880 
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20 
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24 
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The leading Q!^(Za)^-term, given by the B40 coefficient, is 
well known and therefore not included here (for a review see 
e.g. Appendix A of Ref. j^l). The above expression is valid 
for P, D states, and for the normalized difference A„ of 5 
states. The quantity TV is defined in terms of the notation 
adopted in Refs. Il3|| . and the two-loop Bethe logarithm 
6l is defined in Refs. |7, 14]. Although has been deter- 
mined numerically only for 5 states (see Ref. 1 14]), it repre- 
sents a well-defined quantity for all hydrogenic states. The 
logarithmic sum /?4 is given by Eq. 0, with the replacement 
Hji i (T*^ V V . Finally, we have 



^^/?5 = i {a^^V^H-E) In 
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Evaluating the general expression (jSJl for P states, we con- 
firm that BQ2{nP) — A " 7 ■ Furthermore, we obtain the 



results 



. 4 , , - 1 /166 
B,,inP,,2)^-NinP) + -^(- 



4 

B6i{nP3/2) ^-N{nP) + - 



1 



V 405 
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Numerical values for N{nP) can be found in Eq. (17) of 11311 . 
Regarding the nonlogarithmic term Bqq, we fully confirm re- 
sults for the fine-structure difference of P states 1 15]. A fur- 
ther important conclusion to be drawn from Eq. (|8ji is that 
all logarithmic two-loop terms of order {Zoif vanish for 
states with orbital angular momentum I > 2. 

We have also verified that the two-loop result (|8} is consis- 
tent with the normalized 5-state difference A„ for the log- 
arithmic terms Bq2 and Bgi, as derived in Ref. |6] (using 
a completely different method). Evaluating all matrix ele- 
ments in Eq. (|8j, we are now in the position to obtain the 
n-dependence of the nonlogarithmic term, which we write as 
BeoinS)-Beo{lS) = 6L(n5) -6L(15)+yl(n), where ^(n) 



4 



is the additional contribution beyond the n-dependence of the two-loop Bethe logarithm. The result for A{n) is 



'■^8 4 
Mn) - ( i - 3 ln(2) 



[N{nS) - N{1S)] - 



337043 94261 
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45 
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135 n 135 n2 
53 35 419 
60 8^ ^ 120^ 
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129600 21600 n 
53 35 419 
15 2^ 
/37793 16 



ln(2) + + — 
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902609 
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4 16 4 , , 

— + — ) (2) 



C(2) ln(2) 
ln(2) + 8C(2)ln(2) 
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Numerically, A{n) is found to be much smaller than 
&L (nS) — bL{lS), which implies that the main contribution to 
BQo{nS) — Bqo{1S) is exclusively due to the two-loop Bethe 
logarithm. As an example, we consider A{5) = 0.370 042 
and 560(55*) - i?6o(l'5) = 21.2(1.1), where the error is due 
to the numerical uncertainty of the two-loop Bethe logarithm 
6l(5S') (seeRef. Q). 



TABLE II: Theoretical values of the normalized Lamb-shift differ- 
ence A„ = 71"^ AE{nS) — A_E(1S'), based on the results reported 
in this Letter [see Eq. Jl IH . Units are kHz. 



9 
10 
II 



187225 
235070 
254419 
264154 
269738 
273237 
275574 
277212 
278405 
279300 



,70(5) 12 

,90(7) 13 

,32(8) 14 

,03(9) 15 

,49(9) 16 

,83(9) 17 

,90(10) 18 

,89(10) 19 

,21(10) 20 

01(10) 21 



279988.60(10) 
280529.77(10) 
280962.77(10) 
281314.61(10) 
281604.34(11) 
281845.77(11) 
282049.05(11) 
282221.81(11) 
282369.85(11) 
282497.67(11) 



The test of standard model theories and the determination 
of fundamental constants (specifically, of the Rydberg con- 
stant and of the electron mass) provide the main motivations 
for carrying out the QED calculations in ever higher orders of 
approximation. Recently, our knowledge of the ground-state 
Lamb shift has been improved by a fully numerical calcula- 
tion of the two-loop self-energy |8]. However, because of the 
structure of the hydrogen spectrum, the decisive quantity for 
the determination of the Rydberg constant from spectroscopic 
data is the normalized difference A„ of the nS — IS Lamb- 
shift. Elucidating discussions regarding the latter point can be 
found near Eqs. (2) and (3) of Ref. 1 16), and in Appendix A of 
Ref. |3|. Accurate theoretical values for A„ can be inferred 
from the results reported here and are compiled in Table HI] 
The Rydberg constant is currently known to a relative accu- 
racy of 6.6 X 10~^2^ limited essentially by the experimental 
accuracy of the 25* — 8D and 25* — 12D measurements (see 



Table V of |3]). Using the improved theory as presented in 
this Letter, it will become possible to determine the Rydberg 
constant to an accuracy on the level of 10^^^, provided the 
ongoing^experiments concerning the hydrogen 15* — 35* tran- 
sition ilTlilSll reach a sub-kHz level of accuracy. 
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